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Abstract 

We propose the measurement of the decay angular distribution of leptons 
from J/ijj's produced at high transverse momentum balanced by a photon [or 
gluon] in hadronic collisions. The polar and azimuthal angular distribution 
are calculated in the color singlet model (CSM). It is shown that the gen- 
eral structure of the decay lepton distribution is controlled by four invariant 
structure functions, which are functions of the transverse momentum and the 
rapidity of the J/ip. We found that two of these structure functions [the lon- 
gitudinal and transverse interference structure functions] are identical in the 
CSM. We present analytical and numerical results in the Collins-Soper and 
in the Gottfried-Jackson frame. 



1 



The measurement of the angular distribution of lepton's from J/0's provides a detailed 
test of the production and decay mechansim of the cc bound state. So far there has been 
intensive experimental studies of the J/ip production rate and transverse momentum dis- 
tributions in hadronic collisions both at UAl |IJ and CDF ||. However, the observed J/ip 
rate was found to be markedly higher than the predicted one's in ||, where in addition to 
the direct charmonium production also contributions from the x an d the production result- 
ing from B decays were taken into account. It has recently been pointed out in [|J that 
at large transverse momentum of the J/ijj an additional production mechanism comes from 
fragmentation contributions of a gluon or a charm quark into charmonium states. 

In this letter we restrict ourselves to the study of J/i/j produced in association with a 
photon ||. This has several advantages. First of all the experimental signature of a J/ip 
[decaying into an e~e + or p, + pair] and a 7 with balancing transverse momentum is a 
very clean final state. Second, the fragmentation contributions || from radiative \ decays^] 
to J/ip + 7 production are small and the dominant subprocess contributing to a J/ip + 7 final 
state is the gluon gluon fusion process. 

For J/0's produced with transverse momentum p T [balanced by the additional photon] 
one can define an event plane spanned by the beam and the J/ip momentum direction which 
provides a reference plane for a detailed study of angular correlations. The decay lepton 
distribution in J/ip — > in the J/ip rest frame is determined by the polarization of the 
J/i/j. Therefore, the study of the angular distribution can be used as an analyzer of the 
.//^-polarization. It is thus possible to test the underlying J/^-production dynamics [in 
our case the color singlet model (CSM) @] in much more detail than is possible by rate 
measurements alone ||. 

For definiteness we consider the angular distribution of the leptons coming from the 
leptonic decay of J/0's produced with non-zero transverse momentum in association with a 
photon in high energy proton-antiproton collisions: 

p(P 1 ) + p{P 2 ) - J/iP{P) + j(k) +X^T{1) + + j(k) + X , (1) 

where the quantities in the parentheses denote the four-momenta. In leading order of per- 
turbative QCD [O(o^)], J/ip + 7 can only be produced in gg fusion: 

9{pJ + g(Pi) -> mP) + 7(*0 +x. (2) 

Denoting hadron level and parton level quantities by upper and lower case characters, re- 
spectively, the hadron and parton level Mandelstam variables are defined by 

S=(P 1 + P 2 ) 2 , T={P X -P)\ U={P 2 -P)\ (3) 

and 

s = (p 1 +p 2 ) 2 = x x x 2 S , 

t = ( Pl -P) 2 = Xl (T-P 2 ) + P 2 , (4) 
u = (p 2 -P) 2 = x 2 (U-P 2 ) + P 2 , 



* The radiative x,J decays can produce J/ip at both low and high p T , but the photon produced 
will be soft [E ~ 0(400 MeV)]. 
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The angular distribution of the leptons from the J/ip has the general form 



da 3 da u+L 



dpi dy d cos 9 d(p lQn dp 2 dy 



;i + cos 2 #)+ ^A (1-3 cos 2 i 



+ Ax sin 29 cos (ft + —A 2 sin 2 9 cos 20 



(5) 



The angles 9 and in Eq. ([|) are the polar and azimuthal decay angles of the leptons in 
the J/0 rest frame with respect to a coordinate system described below and y (pr) denotes 
the rapidity (transverse momentum) of the J/0 in the laboratory frame. Note that the 
unpolarized differential production cross section denoted by <j u+L is factored out from the 
r.h.s of Eq. (^). The angular coefficients Ai characterize the polarization of the J/0 (see 
below). They are dependent on the choice of the z axis in the J/0 rest frame and are defined 

by the following ratios of helicity cross sections ^da u+L ' L,T ' 1 = ^ 

2 da L 2^2 da 1 4 da T 

A °-da^> M ~d^^ (6) 

The hadronic helicity cross sections dcr ^ I "^ T,/ are obtained by convoluting the partonic 
helicity cross sections sd<J ^ M L ' T ' 7 with the parton densities. One has 

da u +L ,L,T,i sda^ 1 ^ 1 



dpi dy 



J dxidx 2 g hl (x u jj 2 F )g h2 (x 2 , fJ.%) - • ( 7 ) 



Each of the partonic helicity cross sections is calculated in the CSM. The unpolarized dif- 
ferential production cross section is denoted by a u+L whereas <t l ' t,/ characterize the polar- 
ization of the J/0, i.e. the cross section for the longitudinal polarized J/0's is denoted by 
a L , the transverse-longitudinal interference by a 1 , and the transverse interference by a T (all 
with respect to the z-axis of the chosen J/0 rest frame). The results for the helicity cross 
sections a u+L,L,T ' 1 are dependent of the choice of the z axis in the rest frame of the J/0. We 
present explicit results for the Collins-Soper (CS) and the Gottfried- Jackson (GJ) frame: 
In the CS frame |9j the z-axis bisects the angle between P\ and — P 2 

CS: Pi= E x (sin7 C 5,0,cos7c S ) , ^ 
P 2 = E 2 (sin 7 C 5, 0, - cos j C s) , 



with COS7C5 = \j {T - mi )(u-ml) = ]/ 7^ > sin 7cs = -VI -coshes • 

In the GJ frame (also known as t-channel helicity frame) the z-axis is chosen parallel to 
the beam axis 

GJ: P 1 = E 1 (0,0,1), 

P 2 = E 2 (sin7 GJ , 0, cos7 GJ ) , 

with 
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cos 7 G j = 1 - {T _2™){u- m ^) = ll+Zl ' sin 7gj = -V 1 - cos^7 GJ . The beam and target 
energies in the rest frame are E\ = (m 2 — T)/(2m^), E 2 = (m 2 — U)/{2m^) . 

The partonic helicity cross sections in Eq. (|7|) are calculated applying the technique 
described in to the CSM [details of the calculation will be presented elsewhere]. They 
are given by [-R(O) denotes the radial wave function of the bound state]: 



s °cs,gj 167r«a s 2 m^ 2 , , 

-^T= 27s lR{0)l H cs,cjS(s + t + u-m,), (10) 



s 2 t 2 u 2 



with 

cs (t-ml) 2 {u-ml) 2 ^ (u-ml) 2 (s-ml) 2 ^ (s-m 2 J ) 2 (t-m 2 J ) 2, 1 } 

h l ut ( Am V + H u t A (12) 

cs 2(t-ml)(u-m 2 ,){(s-ml) 2 (t-minu-ml) 2+ CS ) ' 1 > 

H CS = 2 H cs > ( 13 ) 

j m^^/sui s(s 2 - ut){t - u) 

CS ~ V2(s-ml) 2 (t -m$) 3 (u-m$ )3 ' t ) 
and for the GJ frame 

H V G r = H u c t\ (15) 



2m 2 ,stu(s 2 + u 2 

iib^j yv — //l^) 4 (-U — uu^ 

rL 



H GJ = 2 H GJ > ( 17 ) 

r m^^/slu(s-u)[s 2 (u-t) + u 2 (s-t)} 
H GJ ~ 7^7" —^77, ___^a,_ 9\2 ■ ^> 



V2(s-m 2 ) 2 (t-m 2 M(u 



m 



As mentioned before, H u+L denotes the matrix element contribution for the production rate 
and was first calculated in the first paper of [0. All other matrix elements correspond to the 
production of polarized J/0's and are given here for the first time. Replacing 16a by 15a s in 
Eq. (|IDD, the results in Eqs. (|10|-|T8|) are also valid for J/0 + g production within the CSM. 



We will now present numerical results for J/ifj + 7 production at the Tevatron collider 
center of mass energy [\^S = 1.8 TeV] including the decay J/ip — > All results are 

obtained using the gluon density parametrization from GRV |TT| with = 200 MeV and 
the one-loop formula for a s with 4 active flavours. If not stated otherwise, we identify the 
renormalization scale n R and the factorization scale \j? F and set them equal to fi 2 = [i 2 F = 
fi 2 R = (rriL +Py(J/V0)- The value for the bound state wave function at the origin |i?(0)| 2 
is determined from the leptonic decay width of J/i/j: T(J/ip —>■ e~e + ) = 4.72 keV, therefore 
|i?(0)| 2 = 0.48 GeV 3 . 
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In Figs. 1 and 2 we show numerical results for the coefficients A% in Eq. ([|) as a function 
of pt{J/iP) and y(Jfif)) in the CS [Figs. l,2(a)] and GJ [Figs. l,2(b)] frame. One observes 
that the coefficients are strongly dependent on Pt{J/4>) both in the CS and GJ frame. As 
mentioned before, the coefficients A Q and A 2 are exactly equal in lowest order in both lepton 
pair rest frames. The angular coefficient A\ is zero in the CS frame for all values of p T (J/ip). 
The reason is that the matrix element for A\ is antisymmetric in u and t [see Eqs. (|14],|6|)1 
and therefore in x\ and x%, whereas the product of the gluon distributions is symmetric 
under the interchange of x\ and x 2 - As a consequence, the rapidity distribution for A\ in 
the CS frame [Fig. 2(a)] is also antisymmetric around y = 0. However, this is different for 
the GJ frame [see Eq. flTB|)]. All coefficients A+ vanish in the limit p (J/ip) — » 0, which can 
be directly seen from our analytical expressions in Eqs. (|TT| - |l8| ,|6|) . 

A similar relation A® Y = A 2 Y was found in LO [0(a s )] in the Drell-Yan process (12 



p + p^V + X— > l + V + X, where V denotes a gauge boson produced at high p T . In JlO 
the complete NLO corrections to the coefficients Af Y are calculated and the corrections 
are found to be fairly small for the ratio's Aj of the corresponding helicity cross section. 
We expect also here, that the LO results for the ratio's A{ in J/ip production are almost 
not affected by higher order QCD corrections. Note also, that the coefficients Ai are not 
dependent on the bound state wave function. The measurement of these coefficients would 
be a sensitive test of the production mechanism of J/0's. 

To summarize: hadronic J/ip + 7 production has been evaluated in the nonrelativistic 
bound state model. In leading order this final state can only be produced by gluon fusion. 
Analytical formula for the decay lepton distributions in terms of four structure functions are 
presented. The angular distribution in the J/ip rest frame is determined by the polarization 
of the J/ip. 
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Figure captions 

Fig. 1 Angular coefficients Aq, A\ and A 2 for J/ip + 7 production and decay in the CS frame 
(a) and GJ frame (b) as a function of the J/ip transverse momentum at y/S = 1.8 TeV. 
No cuts have been applied. 

Fig. 2 Same as Fig. 1 for the y(J/if)) distribution. 
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